Draft version February 1, 2008 

Preprint typeset using LATpX style emulateapj v. 04/03/99 



o 
Q 

(N 



> 

in 
in 

(N 

On 
On 

^ ' 

Or 

6 

c3 



13 



ON THE X-RAY HEATED SKIN OF ACCRETION DISKS 

Sergei Nayakshin 1 

NASA/GSFC, LHEA, Code 661, Greenbelt, MD, 20771 
Draft version February 1, 2008 

ABSTRACT 

We present a simple analytical formula for the Thomson depth of the X-ray heated skin of accretion disks valid 
at any radius and for a broad range of spectral indices of the incident X-rays, accretion rates and black hole masses. 
We expect that this formula may find useful applications in studies of geometry of the inner part of accretion flows 
around compact objects, and in several other astrophysically important problems, such as the recently observed 
X-ray "Baldwin" effect (i.e., monotonic decrease of Fe line's equivalent width with the X-ray luminosity of AGN), 
the problem of missing Lyman edge in AGN, and line and continuum variability studies in accretion disks around 
compact objects. We compute the reflected X-ray spectra for several representative cases and show that for hard 
X-ray spectra and large ionizing fluxes the skin represents a perfect mirror that does not produce any Fe lines 
or absorption features. At the same time, for soft X-ray spectra or small ionizing fluxes, the skin produces very 
strong ionized absorption edge and highly ionized Fe lines that should be observable in the reflected spectra. 

Subject headings: accretion, accretion disks — radiative transfer — line: formation — X-rays: general — 
radiation mechanisms: non-thermal 



1. INTRODUCTION 

X-ray illumination of an accretion disk surface is a problem 
of general astrophysical interest. Since the X-ray heating of 
the disk atmosphere changes energy and ionization balances 
there, the spectra emitted by X-ray illuminated accretion disks 
in any wavelength may be quite different from those resulting 
from non-illuminated disks. It has been known for many years 
(e.g., Basko, Sunyaev, & Titarchuk 1974) that X-ray illumina- 
tion leads to formation of a hot (and often completely ionized) 
X-ray "skin" above the illuminated material. Unfortunately, 
due to numerical difficulties, most previous studies of X-ray il- 
lumination had to rely on a constant density assumption for the 
illuminated gas (e.g., Ross & Fabian 1993; Zycki et al. 1994; 
Matt, Fabian & Ross 1993; 1996; Ross, Fabian & Brandt 1996), 
in which case the completely ionized skin forms only for very 
large ionization parameters (e.g., Ross, Fabian & Young 1999). 

Recently, Nayakshin, Kazanas & Kallman (1999; hereafter 
paper I) have shown that if the assumption of the constant den- 
sity is relaxed, then the temperature and ionization structure of 
the illuminated material is determined by the thermal ionization 
instability, and that the X-ray heated skin always forms on the 
top of the disk (see also Raymond 1993; Ko & Kallman 1994; 
Rozahska & Czerny 1996). Due to the presence of the hot skin, 
the resulting reflected spectra are quite different from those ob- 
tained with the constant density assumption. 

Unfortunately, computations similar to those reported in pa- 
per I are rather numerically involved and time consuming and 
thus are not readily performed. In this Letter, we present an 
approximate expression for the Thomson depth of the hot skin 
which allows one to qualitatively understand effects of the ther- 
mal instability on the reflected spectra. As an example, we ap- 
ply our approximate expression to the observed X-ray Baldwin 
effect (e.g., Nandra et al. 1997b) in the geometry of a cold ac- 
cretion disk illuminated by an X-ray source situated at some 
height above the black hole. 



2. THOMSON DEPTH OF THE SKIN 

We assume that X-rays are incident on the surface of an ac- 
cretion disk whose structure is given by the standard accretion 
disk theory (e.g., SS73). Let us define the hot skin as the region 
of the gas where temperature significantly exceeds the effective 
temperature (see, e.g., Fig. 3b below). As explained in paper I, 
even Fe atoms are completely ionized in the X-ray skin if the 
X-ray flux, P x , is comparable with or larger than the disk intrin- 
sic flux, Fd, and if the illuminating X-rays have relatively hard 
spectra, i.e., photon index T ^ 2. Under those conditions, the 
Compton temperature, T c , is close to ~ few keV. Therefore, we 
can neglect all the atomic processes and only consider Compton 
scattering and bremsstrahlung emission in our analytical study 
of the X-ray skin. 

The Thomson depth of the hot layer, t\, is obtained by inte- 
grating d,T\ = aTn e (z)dz, from z — Zb to infinity, where n e (z) 
is the electron density; z is the vertical coordinate; Zb ~ H 
gives the location of the bottom of the ionized skin; and H is 
the disk scale height. Note that a simpler estimate of t\ can be 
obtained if one assumes that the gas temperature is equal to the 
Compton one, and that the density law follows a Gaussian law 
(see Kallman & White 1989). 

In our calculations of t\, we will neglect by the reprocessing 
of the radiation as it penetrates through the skin. This limits the 
applicability of our results to t\ 1 if the angle, 9, that the X- 
ray radiation makes with the normal to the disk, is not too large, 
i.e., that /z.j = cos 6 is not too small. In the case when fa -C 1, 
the spectral reprocessing cannot be neglected for t\ £ fa. For 
simplicity, we will assume \ii £ 1 and postpone the treatment 
of large incident angles to a future publication. 

If -Pbot is the gas pressure at z = Zb, then we can define di- 
mensionless gas pressure as P» = P/Pbat, where P — p/pkT 
and p is the gas density, and also dimensionless gas tempera- 
ture by = T/T c (T c is the Compton temperature, see be- 
low). The Compton-heated skin then has P* < 1, whereas the 
colder material below the skin has P* > 1. As shown in paper 
I, the incident X-rays do not affect the hydrostatic balance in 
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the ionized skin because the main source of opacity is Comp- 
ton scattering (see Fig. 5b in paper I and also Sincell & Rrolik 
1997). With this, one can re-write hydrostatic balance equation 
in terms of dimensionless variables as 



<9P* 
dx 



,-P* £-C 
J T* A 2 



(1) 



where x = z/H, £ is the ratio of the disk midplane radia- 
tion pressure to the total pressure £ = Prad(0)/(P ra d(0) + 
Pg as (0)), and A is the scale-height of the skin in units of H: 
A = (2kT c R 3 /GM/iH 2 ). Expressing the electron density as 
n e = p/n e , one can arrive at 



Tl 



A 4 c -fbot lx 



/'-<:: 



P* 
TV 



dy = t W(y b ) , (2) 



where l x is the compactness parameter of the illuminating X- 
rays (see, e.g., equation 20 in paper I); 9 C = kT c /m e c 2 ; we 
also defined y = (x — £) /A and y b — (zb/H — QJ A; finally, To 
is the expression before the sign of the integral. The integral in 
the equation (Q) is designated as W(yb) and often (see below) 
turns out to be of order unity. Krolik, McKee & Tarter (1981, 
§IV2b; KMT hereafter) showed that for a completely ionized 
gas the energy balance equation can be written as 



T 



(3) 



where S* = S/S^ c is the pressure ionization parameter normal- 
ized by the "inverse Compton" ionization parameter S,- c , which 
is given by equation (4.5) of KMT, re-written in order to take 
into account the difference in definitions of 3 = F x /cP (the 
one used here) and 3 = Fi on /riHkT (KMT; Fi on is the X-ray 

flux between 1 and 10 3 Ry): 5 ic = 0.47 * T~ 3/2 ; T 8 here is the 
Compton temperature in units of 10 8 Kelvin. 

The solutions presented in paper I possess the property that 
the transition from the Compton-heated to cooler layers occurs 
at the upper bend of the S-curve (e.g., point (c) in Fig. 1 of 
paper I). In that point, dT/dE = oo. Using this condition, 
one can show that the transition happens at T* = 1/3 where 

„ o3/2 

s = ^i C . Therefore, 



Pbot = 



S,^ 1 — = 3.2 x 10~ 2 Tf /2 — 
3\/3 c 1 c 

9I6 

G(r)^ Z —[l-{i/rfl 2 fr-\ 



(4) 
(5) 



where Ti = kT c /l keV; r = R/R s , R s = 2GM/c 2 is the 
Schwarzschild radius; m is the dimensionless accretion rate 
(m = 1 corresponds to Eddington luminosity for the accretion 
disk) and < / < 1 is the coronal dissipation parameter (e.g., 
Svensson & Zdziarski 1994). The maximum of function G(r) 
occurs at r = 16/3 where it is equal to unity. We also approx- 
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/i = to p /2, and W(y b ) ~ W(0) = 1.23 for 



the following reason. If the vertical pressure profile is given by 
P(z) = P(0) exp[— (z/H) 2 ] for the gas-dominated disk, the 
location of the temperature discontinuity is 



= H In 1 / 2 



P(0) 



P 



bol 



(6) 



The function ln 1 ^ 2 (i) is a very slow, monotonically increas- 
ing one: for example, its value is 1.52 for t = 10 and 3.72 
for t = 10 6 . Therefore, for most realistic situations, z b is not 
very much larger than H. As is easy to check, A ^> 1 for 
gas-dominated disks, and thus y b -C 1. For radiation-pressure 
dominated disks, z b is nearly equal to H (see §3.4 in paper I), 
and hence y b = (1 — C)/A <SC 1 in this case as well. To summa- 
rize this statement in words, the vertical extent of the ionized 
skin is always large enough that the exact location of the inner 
boundary is unimportant, for either gas or radiation-dominated 
disks. 

A cautionary note is in place. So far we assumed that the 
spectrum of the ionizing radiation does not change with the 
depth into the illuminated material. This is true only for op- 
tically thin situations, i.e., when t\ ^ 1, but for larger optical 
depths the Compton temperature is smaller than the one corre- 
sponding to Ti = 0. The decrease in the local value of T\ dic- 
tates a corresponding decrease in the value of Pbot according 
to equation (0), and therefore the transition to the cold equilib- 
rium solution will always happen at a smaller t\ than equation 
(||) predicts. To reflect this fact, Kallman & White (1989) intro- 
duced TTh crit ~ 1 — 10, defined to be the maximum of t% . Note 
that the behavior of Th with Compton temperature, X-ray flux 
and radius given by equation (Q) coincides with that obtained 
by these authors (Tics on the second line of their equation [1] 
should be in power + 1 rather than — 1 /2 [T. Kallman 1 999, pri- 
vate communication]). 

With the understanding that for t\ it 1 our analysis may not 
be expected to be very accurate, we choose to fix rxh crit at 
a value of 3. This number is motivated only by the conve- 
nience - in paper I, due to the large volume of calculations, 
we have limited the Thomson depth of the illuminated gas that 
can be strongly ionized to tt = 3. This choice did not lead 
to any spurious results because the reprocessing features (from 
the cold material below the skin) become negligible already at 
tt = 2 (see paper I). In this paper, we use the following simple 
modification to the Thomson depth of the hot skin: 



Th 



i + n/3 



(7) 



We can now compare our results with values of Th numer- 
ically calculated in paper I. There, we conducted two types 
of tests. For radiation-dominated disks, we used fully self- 
consistent formalism and thus we can directly compare results 
of §5 in paper I to results of equations (Q) and ([7J). For pur- 
poses of isolating dependence of the reflected spectra on one 
parameter at a time and yet covering a large parameter space, 
we also conducted several tests where we artificially varied the 
"gravity parameter A" (see §4 in paper I). In order to com- 
pare results of these latter tests, one should use the same ap- 
proach that led us to equation (|j), and also use the same gravity 
law as used in paper I with (note that in our current formula- 
tion, A = 46» C J- 1 A- 2 ). This yields ti = 0.73 Ti (Ix/A) 1 / 2 . 
Further, since A is now expressed as A 2 = 46> c v4~ 1 1" 1 , the 
value of yb is not necessarily small, and thus W(y b ) needs 
to be evaluated exactly. Finally, to locate the position of the 
lower boundary of the illuminated layer (which occurs almost 
at the same height as the temperature discontinuity - see Fig. 
5d in paper I for an example), we must substitute Pbot hi 
equation (^) on the value of gas pressure at Thomson depth 
Tmax = 4, as used in paper I, which is approximately equal to 
r max AF x /c = 4AF x /c. 
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The results of such a comparison are summarized in Figure 
(1), where we show Th obtained in this paper versus that ob- 
tained numerically in paper I. The value of the Compton tem- 
perature (i.e., Xi in eq. ||) is taken to be the temperature of the 
first zone in the temperature profiles shown in the correspond- 
ing Figures in paper I. We note that the deviation of our approx- 
imate expression from "exact" results is less than ~ 20% even 
though we covered a wide range of physical conditions, i.e., 
strong to weak illumination limits, different indices of the inci- 
dent X-ray radiation, and the disk itself is either gas or radiation 
dominated. 

3. "LAMP POST" MODEL 

As an application of our methods, we choose to analyze the 
model where the X-ray source is located at some height h x 
above the black hole (the "lamp post model" hereafter). Iwa- 
sawa et al. (1996) reported observations of iron line variability 
in Seyfert Galaxy MCG-6-30-15, and pointed out several prob- 
lems connected with theoretical interpretation of these observa- 
tions. Reynolds & Begelman (1999) showed that the accretion 
flow within the innermost stable radius may be optically thick 
and thus produce fluorescent iron line emission in addition to 
such emission from the disk itself. These authors argued that 
the line emission from within the innermost stable radius may 
be important for the interpretations of the observations of Iwa- 
sawa et al. (see also Dabrowski et al. 1997). Reynolds et al. 
(1999) studied response of the iron line profiles to changes in 
the X-ray flux (iron line reverberation). In this paper we will 
not discuss the region within the innermost stable orbit for the 
reason that the properties of the accretion flow there are not 
well constrained. In particular, the hydrostatic and energy bal- 
ance equations do not necessarily apply since the hydrostatic 
and thermal time scales can be longer than the in-fall time. 

We will assume a non-rotating black hole and that all the X- 
rays are produced within the central source. We also neglect all 
relativistic effects in the present treatment. The projected X-ray 
flux impinging on the disk is 



L x hz 



Att(R 2 +h? x fl 2 



(8) 



Let us define r\ x as the ratio of the total X-ray luminosity to 
the integrated disk luminosity of the source, i.e., rj x = L x jL^. 
Using equation (||), one obtains 



n ~ 27.2 ^ r) x Ti [1 + (h x /Rf]- z l 2 r- 3 / 2 tin (9) 



(one still has to use equation [|7j] to get the final answer for Th). 
For illustration, we choose a value of h x — 6R S - For T — 1.9, 
a typical value for Seyfert Galaxies, the Compton temperature 
kT x ~ 7.6 keV. As discussed in paper I (§4.6), the Compton 
temperature at the surface of the disk depends on the cosine 
of the X-ray incidence angle, /!;, and the ratio F^/F^ approxi- 
mately as: 



T ~ T 

1 r — J- cr. 



Hi a/3 — 



Fa 



(10) 



Figure (2) shows the Thomson depth of the skin as a function 
of radius for several values of rh. The parameters in Figure (2) 
are chosen to be: (a) T = 1.9, r) x = 1 ("X-ray strong" case); 
(b) T = 1.9, r] x = 0.1 ("X-ray weak" case), and (c) V = 2.3, 



r\ x = 1. Figure (3) shows the angle-averaged reflected spectra 
and temperature profiles of the hot layer for the three cases just 
considered and the accretion rate equal to the Eddington value 
(rh = 1) for r = 10. The details of numerical methods with 
which these spectra were obtained are described in paper I. 

If Compton temperature T c ^ 1, the illuminated gas is nearly 
completely ionized, and the line is emitted almost exclusively 
from the material below the skin. For this reason, it turns out 
that local emissivity of the line is negligible when Th >, 1 (see 
paper I). Thus, the strength of the iron line will be decreasing 
with rh for the X-ray strong case (Fig. 2a) when L x £ few 
percent of the Eddington value. The skin is thickest for smaller 
radii, and therefore the broad iron line component will decrease 
first. The narrow line component (emitted farther away from 
the black hole) will also decrease with X-ray luminosity, but 
considerably slower than the broad line. In addition, if some 
of the line comes from a putative distant obscuring torus (Kro- 
lik, Madau & Zycki 1994), or the disk has a concave geome- 
try (Blackman 1999), then the narrow component will decrease 
even slower. Thus, our qualitative description of the behavior 
of the iron line EW with luminosity for X-ray strong sources is 
consistent with Figure (3) of Nandra et al. (1997b), suggesting 
an explanation to the X-ray Baldwin effect. 

The X-ray weak case (Figs. 2b & 3b) is different in two re- 
spects. Firstly, the Thomson depth of the hot skin is smaller 
at a given accretion rate compared with the X-ray strong case. 
Most importantly, however, the skin is not "that hot". Namely, 
the skin temperature is only ~ 0.3 keV in the inner disk, and 
some of the iron ions are not completely ionized. For that rea- 
son, it turns out that the iron line centroid energy is close to 6.7 
keV, and a very deep absorption edge appears. This edge is in 
fact much stronger than the one resulting from a neutral mate- 
rial. Hence it is possible that this relatively cold X-ray heated 
skin can be unambiguously detected in spectra of real AGN. 

Similarly, the soft incident X-ray spectrum leads to a rela- 
tively cool skin because T x ~ 1 keV only (see Fig. 2c & 3c). 
As in the case r\ x = 0.1, a strong absorption edge is observed. 
In addition, the 6.7 keV Fe line is stronger, with EW of 65 eV. 
Therefore, such skin may also be detectable if it exists. 

Comparing the value of Th resulting from analytical formu- 
lae with those seen in Fig. (3b), one notes that deviations are 
as large as ~ 50%. These relatively large deviations exemplify 
the fact that our equations (Q) and (Q) are good approximations 
only to cases with strong X-ray flux and hard incident spectra 
(i.e., F x >, Fd and V < 2). When one of these two conditions 
is not satisfied (and it is the case for all three cases presented 
in Figs. 2 & 3) the Compton heated layer is cool enough for 
atomic processes to provide additional sources of heating and 
cooling beyond Compton and bremsstrahlung processes, so that 
the energy equation ( |) is not obeyed. However, our results can 
still be used as an order of magnitude estimate of Th in the latter 
cases. 



4. SUMMARY 

In this Letter, we derived an approximate expression for the 
Thomson depth of the hot completely ionized skin on the top 
of an accretion disk illuminated by X-rays. Our results are only 
weakly dependent on the a priori unknown a-viscosity param- 
eter (because it enters through the boundary conditions). This 
allows us to reduce the uncertainty in documenting predictions 
of different accretion theories with respect to the iron line pro- 
files and the strength of the X-ray reflection hump. Using the 



4 



"lamp post model" geometry as an example, we showed that 
an inner part of an accretion disk may have a Thomson-thick 
ionized skin. Under certain conditions (F x £ Fd and T £ 2), 
the X-ray heated skin may act as a perfect mirror for photons 
with energies below ~ 30 ke V. The physical cause of this is that 
Compton scattering in the skin prevents X-rays from penetrat- 
ing to the deep cold layers that are capable of producing fluo- 
rescent line emission as well as setting other marks of atomic 
physics. 

We note that because the ionized skin is thickest in the in- 
ner disk, the observed absence or deficit of the relativistically 
broadened line and other reprocessing features in some systems 
(e.g., Zycki , Done & Smith 1997, 1998), which was inter- 
preted as a possible evidence for a disruption of the cold disk for 
small radii, may also mean that the "cold" disk is still present 
up to the innermost radius, but the skin effectively shields it 
from the X-ray illuminating flux. It is interesting, however, 
that the presence of the skin becomes apparent in systems that 
have f x ^ F,i or T ^, 2, because the Fe atoms are not com- 
pletely stripped of their electrons and thus produce strong ion- 



ized edges and lines. We believe these predictions should be 
testable observationally with current X-ray missions such as 
Chandra, Astro-E and XMM. Also note that for a patchy corona 
model of accretion disks (e.g., Haardt, Maraschi & Ghisellini 
1994), the Thomson depth of the hot skin is always larger than 
the one found here, since the ratio F x / Fj is larger. Finally, 
the presence of the ionized skin is important not only for the 
X-rays, but for other wavelengths as well (e.g., in studies of 
Lyman edge of accretion disks, correlation of optical/UV light 
curves with X-rays). 

A shortcoming of this paper is that we neglected the changes 
in the ionizing continuum due to scattering in the skin, which, 
strictly speaking, limits the applicability of our results to opti- 
cally thin situations ,$ 1 and incident angles not too far from 
normal 2 . We plan to address this issue in future. 

The author acknowledges support from NAS/NRC Associ- 
ateship and many useful discussions with D. Kazanas, T. Kail- 
man and M. Bautista. The anonymous referee is thanked for 
several constructive suggestions. 
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FIG. 1. — Comparison of the analytical approximation to the Thomson depth of the hot skin and the numerical results presented in paper I. Filled circles show 
tests with r = 1.9, but with varying gravity parameter A (§§4.2-4.4 in paper I); crosses show tests with A = 0.3 and with T varying between 1.5 and 2.3 (§4.7 
in paper I); and stars show the radiation-dominated case (described in §5 - dotted and long-dashed curves in Fig. 13 in paper I). The solid curve shows the line 
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FIG. 2. — Thomson depth (solid curves) of the hot skin as a function of radius for different accretion rates (whose values are shown next to corresponding curves). 
The Compton temperature is shown by the dotted curve. See text for additional details. 
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Fig. 3. — Reflected spectra (a) and temperature profile (b) computed for the lamp post geometry at r = 10 with accretion rate fa = 1. The dashed, dotted and 
solid curves correspond to models shown in Fig. (2) [a], [b] and [c], correspondingly. The solid curve in [a] was scaled down by factor of 2 for clarity. 



